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Abstract 

Magnetic reconnection (MR) in the Earth's magnetotail is usually followed by a system-wide 
redistribution of explosively released kinetic and thermal energy. Recently, multi-spacecraft ob- 
servations from the THEMIS mission were used to study localized explosions associated with MR 
in the magnetotail so as to understand subsequent Earthward propagation of MR outbursts dur- 
ing substorms. Here we investigate plasma and magnetic field fluctuations/structures associated 
with MR exhaust and ion- ion kink mode instability during a well-documented THEMIS MR event. 
Generation, evolution and fading of kink-like oscillations is followed over a distance of ~ 70,000 km 
from the reconnection site in the mid- magnetotail to the more dipolar region near the Earth. We 
have found that the kink oscillations driven by different ion populations within the outflow region 
can be at least 25,000 km from the reconnection site. 
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Collisionless magnetic reconnection (MR) necessitates formation of a thin current sheet 



in the magnetotail 



]J. The current sheet can develop instabilities, including the kink mode. 



Kink-like flapping fluctuations were extensively observed by the four Cluster spacecraft 

aa. 

Accordingly, the internally driven kink-like waves propagate from the center to the flanks 
of the magnetotail with a typical period of a few tens of s. Kink motions can be induced, 
for example, by the magnetic double gradient mechanism or by ion bulk speed shears, 
even under bifurcated current sheet configurations jlj]. The instabilities can also exhibit a 
large growth rate driven by relative streaming between separate ion populations. Though 
full particle and hybrid simulations with open boundary conditions 6] revealed many salient 
features and observable signatures of the ion-ion kink mode driving and evolution of asso- 
ciated fluctuations are not fully understood. Recent full particle simulations and analysis 
indicate that the ion temperature in the reconnection outflow is proportional to proton bulk 
speed and ion species mass These results allow us in this paper to distinguish the 
separate ion populations essential for driving the ion-ion kink. Moreover, we analyze the 
evolution of MR outflow and kink-associated fluctuations using unique observations from 
the THEMIS mission. The THEMIS mission, with its fleet of five identical spacecraft, was 
designed to study substorms during which a system-wide reconfiguration of the magneto- 
sphere occurs. THEMIS timing analysis (comparison of particle and field MR signatures at 
the spacecraft positions) shows that system-wide substorm activity is initiated locally and 
the likely triggering mechanism is MR in the mid- magnetotail {§]]. 

We analyzed the field and plasma fluctuations associated with MR from 4.8 to 5.3 UT on 
26 February 2008 (hereafter, decimal hours will be used). Substorm- related timing analysis 
of this event, together with ground-based (geomagnetic observatory), auroral and magne- 
totail signatures were described in great detail in [8|. We used magnetic data with 3s time 
resolution from the THEMIS FGM experiment 9| and 3s resolution ion moments computed 
from both ESA (energy range: eV - 25 keV) 10] and SST (energy range: 25 keV - 6 MeV) 
instruments. The THA spacecraft (TH - THEMIS; A-E, specific spacecraft), situated at 
~5.5 Re (Earth radii), did not observe significant bulk flows. Because the observations 
from THD and THE, which were separated by ~ 0.8 Re, were similar, data from THA 
and THE will not be used further. The probes THB, THC and THD were aligned along 
X (from -21 to -11 Re-, see positions indicated in Figure 1) within 3.9 < Y < 4.6 Re and 
-3 < Z < -2 Re in the GSM system of coordinates. Three probes were lined up from the 
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mid-magnetotail to near-Earth space. Since the average Y-Z cross-section of MR-driven 



bulk flows is 1-3 Re 



the same Earthward flows could be observed by the spacecraft. 



Figure 1 shows the magnetic field components Bx, Bz and the bulk speed component Vx in 
three vertically aligned subplots for each spacecraft. The subplots in the third row indicate 
that THB saw MR-associated flow reversal (before and after 05.0 h UT: tailward - Earth- 
ward flows, — Wx), whereas the other spacecraft observed predominantly Earthward flows. 
The flows reached maximum values of \Vx\ > 500 km/s. The flow reversal is interpreted as 
movement of the reconnection X line through THB's position. Because THC approached 
the plasma sheet from the lobe (Bx < —20 nT), it did not observe the flow onset. The 
Earthward flows were associated with fluctuations in the Bx and Bz magnetic components 
(top two rows in Figure 1). Although the largest amplitude fluctuations in Bx occurred 
closer to the MR site (at THB and THC), they remained significant at THD. Since Bx was 
frequently changing sign, the observed fluctuations could at least partially be attributed to 
a flapping current sheet. The increase in fluctuations and magnitude of Bz (second row 
in Figure 1) during Earthward flows can be attributed to flow-driven dipolarizations of the 
magnetic field that are more significant closer to the Earth (at THD, last subplot). 

Now we examine the orientation of current sheet oscillations associated with the Bx 
sign change. The current sheet orientation is estimated using minimum variance analysis 
(MVA), which is suitable for single spacecraft measurements. Assuming a 1-D boundary, the 
eigenvector (ni) corresponding to the smallest eigenvalue (A;) computed from the covariance 



matrix of magnetic field is taken as the boundary normal 12[ . At THB and THC the flow 



associated current sheet normal vectors pointed predominantly in the Y-Z directions (n3y vs. 
n3z pairs lay on a circle, not shown), n3x was close to zero. This behavior is characteristic of 
the Y-Z kink-like mode under tilted current sheet conditions [ijj . Nevertheless, this finding 
does not exclude the occurrence of a mixture of kink- and sausage-mode fluctuations. 

MVA normals showed that at THD, the Y-Z kink signatures were largely lost, possibly 
because of different types of interactions during dipolarization and flow braking closer to 
the Earth. 

Let us now investigate high-speed flow intervals more closely. The first column in Figure 2 
shows superimposed Vx and Bx fluctuations for each spacecraft. The axis values for Bx are 
the same for the first and second columns in the figure. The second column shows Bx as a 
function of plasma (3. The horizontal dashed lines represent Bx = nT; the horizontal thick 



3 



[nT] 



10 


-10 
-20 
20 
10 


-10 



[nT] 



V 



x 



500 



[km/s] u 

-500 



February 26, 2008 























; Jli 

UT 


[h] 


ill 



4.8 4.9 5.0 5.1 5.2 4.8 4.9 5.0 5.1 5.2 4.8 4.9 5.0 5.1 5.2 5.3 

THB (-21 R_) 



THC (-17 R E ) 



THD (-11 R E ) 



FIG. 1: Bx, Bz magnetic and Vx bulk speed GSM components observed by THEMIS: THB 
(left), THC (middle), THD (right). 

lines, high-speed flow intervals. The correlation coefficients between Bx and Vx fluctuations 
were 0.66 and 0.68 at THB and THC, respectively. The correlations between current sheet 
flapping motions (±Bx) and bulk speed (Vx) fluctuations can be partially explained through 
a spatial effect: plasma (3 increases towards the center of the current sheet and towards the 
central part of the plasma flows. In fact, the condition (3 > 2 indicates that a spacecraft 
is merged to a high-speed bulk flow 11] . As the probe moves through the current sheet, it 
can intersect different parts of the flow again and again. The correlation between Bx and 
Tj (not shown) emerges for the same reason. The right column shows superimposed Vx and 
Ti/10 fluctuations for each spacecraft. The correlation coefficients are higher between the 
plasma parameters 0.91 and 0.86 than between the magnetic field and plasma parameters 
0.66 and 0.68 for THB and THC, respectively. The fluctuations at THD are uncorrelated, 
since the spacecraft are in the center of the plasma sheet in the dipolarization region, where 
plasma (3 is changing less. Along with MVA analysis results, the figure shows that magnetic 
and plasma fluctuations evidencing of kink oscillations span a distance of at least ~ 4i?£ 
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FIG. 2: (Color). Correlations between Vx and Bx (left column); dependence of plasma f3 on Bx 
(middle column); correlations between Vx and Tj/10; THB (top); THC (middle); THD (bottom). 

(from THB to THC), but do not reach the position of THD. A similar large-scale extension 
of kink-like oscillations in the magnetotail was observed only once during a conjunction of 
Cluster and Double Star satellites. \}\ . 

An important question is what drives speed and magnetic kink-like fluctuations during the 
high-speed MR outflows at THB and THC. Global MHD simulations show that directional 
changes in solar wind speed (mainly in Vz) can induce neutral sheet flapping in the near 
and middle magnetotail with a delay of 10 to 15 minutes [ljj] We have checked time-shifted 
(to the distance of the Earth) ACE solar wind data between 04.66 and 05.66 UT and found 
the following mean and standard deviations: Vx = —385 ± 2 km/s, Vy = —10 ± 3 km/s, 
V z = 32 ± 3 km/s, B x = 0.3 ± 1 nT, B Y = 2.2 ± 1 nT and B z = 1.2 ± 0.3 nT. Since the 
changes are insignificant, flapping current sheet motions are unlikely to be driven by the solar 
wind. Moreover, the current sheet is quiet between tailward and Earthward flows ( top left 
subplot in Figure 2). Intense current sheet oscillations and fast flows arose simultaneously 
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FIG. 3: (Color). Comparison of the time evolution of plasma residuals and the flapping magnetic 
field at THB (top) and at THC (bottom). 

If the solar wind were driving the current sheet, strong Bx oscillations would be present 
between the tailward and Earthward flow regions, as well. Current sheet oscillations due to 
finite By in the neutral sheet o r a number of inabilities Q ean oecur during the entire 
reconnection interval [l?]]. We are not aware of any experimental work that clearly identifies 
a single mechanism driving current sheet oscillations in space. We show here that different 
ion populations exist near MR (at THB), and strong interactions between the magnetic 
field and plasma outflow are present near and at least 4 Re away from the MR site. These 
interactions could be responsible for current sheet kink oscillations during the fast flows. 
Since the correlations between Bx-, Vx and during fast flows are lost at THD (see the 
bottom row of subplots in Figure 2), we consider only the fluctuations at THB and THC. To 
remove the spatial effects due to multiple crossings of the current sheet and plasma flows (for 
which the changing plasma (3 serves as a proxy measure), we evaluated the linear fit between 
Bx, Vx and B x , T% during fast flows in a least-squares sense, respectively. We suppose, the 
residuals of the fit correspond to the fluctuations of plasma parameters not influenced by the 
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FIG. 4: (Color). Cross- wavelet spectra between Vx res and Bx at THB (top) and THC (bottom); 
colors indicate cross-wavelet power; thick contours represent 5% significant level against red noise; 
the arrows correspond to the cross- wavelet phase: with in-phase (anti-phase) pointing right (left), 
and Vx res (Bx) leading Bx (Vx res ) by 90° pointing down (up). The lighter shade indicates regions 
where edge effects might distort the spectra. 

flapping current sheet. The physical significance of residuals is also supported by the higher 
correlation between plasma parameters than between plasma and Bx fluctuations; that is, 
plasma parameter fluctuations are not determined exclusively by flapping. Figure 3 shows 
the relationship between the residuals from the linear fit (Vx res , Ti res ) and Bx at THB (top 
panels) and THC (bottom panels). On the right, Bx, Vx res and Ti res /10 time series are 
depicted. On the left, the 3D panel shows the the squares of the speed and temperature 
residuals as functions of time. 

There are several position-dependent items of interest. Close to MR (at THB) the residu- 
als (top-left in Figure 3) show two well-correlated plasma populations, I and II. Population I, 
characterized by its low speed and high temperature, appeared before ~ 05.09 UT; Pop- 
ulation II, characterized by its high speed and low temperature, appeared after ~ 05.09 
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UT. Recent full-particle simulations indicate that ions crossing the MR outflow boundary 
can behave like non-adiabatic pickup particles, gaining an effective thermal speed depend- 
ing on the bulk speed of the reconnection exhaust. The total ion temperature depends on 
the square of the outflow speed and the mass of ion species Q]. Cluster analysis has al- 
ready revealed that ion motion is usually non-adiabatic in the flapping current sheets fis| . 
This particle pickup mechanism represents a possible explanation of the observed speed- 
temperature correlation. The combination of Y-Z like kink motions (see above) and the 
deduced change in ion populations indicate the occurrence of ion-ion kink instability. More- 
over, we show that the dynamics of associated current sheet flapping changes when the 
plasma parameters change from population I to population II. The correlation coefficient 
between the residuals is r(Vx res , Ti res ) = 0.87, and the effect of flapping is largely removed, 
since r(Vx res , Bx) ~ r(Ti res , Bx) ~ 0.02. Nevertheless, the cross-wavelet power between 
Vx res (or Ti res , not shown) and B x is significant (top subplot in Figure 4). The cross- 
wavelet plot associated with population I plasma (before ~ 05.09 UT) displays significant 
power over a range of periods (~ 4-20 s), while longer periods (~ 15-60 s) are associated 
with the population II plasma (after ~ 05.09 UT). An inspection of Bx data (right in Figure 
3) also indicates that the character of Bx fluctuations changes when the border (a vertical 
line at 05.09 UT) between plasma populations I and II is crossed. Shorter-period magnetic 
fluctuations (~ 10-70 s) are replaced by longer-period (~ 70-120 s) ones. The longer periods 
represent typical periodicities found in 3D full particle and hybrid simulations during the late 
stage of ion-ion kink instability 6j and a characteristic period observed by Cluster during 
current sheet kink oscillations Il9j. Simulation results also show that short period ion- ion 

n 

kink oscillations occur first, evolving later into longer period modes [6J. The transition from 
shorter- to longer-period fluctuations is visible in both Bx (Figure3) and cross-wavelet spec- 
tra (Figure 4) at THB. At THB and THC positions, the total magnetic field (not shown) 
strongly fluctuates between 2-17 nT. The corresponding proton gyroperiod varies between 
4-20 s, overlapping the range of observed shorter-period fluctuations. 

At the position of THC the correlation between the plasma parameter residuals and mag- 
netic field (r(Vx res , Bx) ~ r(Ti res , Bx) ~ 0.04) is negligible. At a distance of 4 Re Earth- 
ward from MR, the interconnection between residuals is weaker (r(Vx res ,Ti res ) = 0.68). 
Though THC is in the flow for less time, no differing plasma populations are visible, and the 
cross- wavelet Vx res , Bx plot shows common power at discrete periods of ~8 and 25 s. Bx 
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shows the largest power over periods of 100-200 s (not shown). The decreasing correlations 
and increasing periods at THC signal that the ion-ion kink reached its saturation or even fad- 
ing phase. The wavelet phase, which would reveal the in-phase or out-of-phase features [20I 
associated with the quantities in the cross-wavelet plane, exhibits knotty, period-dependent 
patterns at THB (the direction of arrows in Figure 4). The pattern is less complex at THC. 
Nevertheless, it seems to be impossible to identify a clear wavelet phase between plasma 
and magnetic field fluctuations. We can only say that the cross-wavelet fluctuations exhibit 
compound multi-scale interactions and the cross-power of these interactions is changing with 
time and position. 

In summary, we have identified reconnection outflow-associated plasma and magnetic 
field correlations in the presence of ion- ion kink mode oscillations of the current sheet. The 
evolution and fading of the kink mode with distance are connected with decreasing speed- 
temperature correlations, therefore with decreasing efficiency of pickup ion feeding. These 
findings suggest that the simultaneous occurrence of multi-scale interactions between the 
flapping magnetic field and the plasma flow can partially be responsible for the gradual 
disappearance of kink mode oscillations. Further analysis is needed, however, to establish 
causal connections between driving and dissipation mechanisms for the observed correlated 
fluctuations. 
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